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Relying only on first principles, we derive a master equation of Lindblad form generally applica- 
ble for adiabatically time dependent systems. Our analysis shows that the much debated secular 
approximation can be valid for slowly time dependent Hamiltonians when performed in an appro- 
priate basis. We apply our approach to the well known Landau-Zener problem where we find that 
adiabaticity is improved by coupling to a low temperature environment. 
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The adiabatic theorem [TJ [2] is an approximation de- 
rived in the very early days of quantum mechanics. 
Its application requires that the time dependence of an 
Hamiltonian H{t) is slow compared to the internal time 
scales, as quantified by a small adiabatic parameter A. 
The adiabatic theorem then states that the time depen- 
dence of the Hamiltonian only induces transitions be- 
tween different instantaneous eigenstates of the order of 
A 2 . For an analytic Hamiltonian, this result can be im- 
proved [31 E] to transitions of order of e -1 /" 4 . 

The adiabatic theorem is of great significance in many 
traditional fields of quantum mechanics, like molecular 
physics, where it is a prerequisite for the use of the 
Born-Oppenheimer approximation [5], and quantum op- 
tics where it is used, e.g. for stimulated Raman adiabatic 
passage [5] and the related electromagnetically induced 
transparency [7] to slow down light to near standstill. 

Renewed interest in the theorem came with Berry's 
discovery of geometric phases in cyclic adiabatic sys- 
tems [8 ( 9J. Berry's phase soon got generalized to non- 
commuting geometric rotations within a degenerate sub- 
space, so called holonomic transformations [TOj. As these 
turned out to be robust to some errors in the control of 
experimental parameters, holonomies received a lot of at- 
tention in the context of quantum computing [111 112) . 

An entirely different approach to quantum computing, 
called adiabatic quantum computing [T5J [TJ] , also heav- 
ily relies on the adiabatic theorem. The aim here is to 
find the ground state of a complicated Hamiltonian i?fi na i 
of interacting qubits. This is achieved by first preparing 
the system in the ground state of an easy Hamiltonian 
initial (e.g. a non-interacting Hamiltonian). The Hamil- 
tonian is then adiabatically changed to i?fi ri ai and the 
adiabatic theorem guarantees that the final state is the 
desired ground state. This process can be used for quan- 
tum computing because it was shown that knowing the 
ground state of an appropriate Hamiltonian is equivalent 
to knowing the solution of the desired calculation. 

The minimal time required for an adiabatic transfor- 
mation is limited by A <C 1 and about an order of mag- 
nitude larger than typical non-adiabatic transformations. 
Even weak coupling to an environment might therefore 
significantly influence the dynamics of the system. In 



particular in respect to quantum computing, where high 
fidelity transformations are needed, decoherence is a ma- 
jor obstacle. It is of great importance to understand the 
interplay of decoherence and adiabaticity. 

Adiabatic systems weakly coupled to an environment 
were already studied in the literature in connection 
to Berry's phase |15H19| . adiabatic quantum comput- 
ing [201 HI], molecular magnets [25], biological sys- 
tems [23], as well as Landau-Zener transitions [24]. In 
some of these the authors postulated a Markovian master 
equation (ME) of Lindblad form similar to time indepen- 
dent systems. However, it is well known that for time 
dependent Hamiltonians the secular approximation can 
cause problems and hence a ME of Lindblad form has 
to be used with care. While f22j [23] derived a Lindblad 
ME from first principles, they neglected terms of order A 
in the dissipative part of the ME. Although reasonable 
for the applications these authors had in mind, in the 
context of Cooper pair pumping it was shown that terms 
of order A have to be included to find the correct trans- 
ferred charge [25H28] . In general, because the adiabatic 
theorem is correct up to e~ x l A it is desirable to derive a 
ME to the same order. This is done below. 

The description of quantum systems by Lindblad MEs 
is standard for time independent systems 29J and is de- 
sirable for two main reasons: First, quantum trajectory 
theory |29j provides not only a neat physical intuition 
for dissipative processes, but it is also an efficient numer- 
ical method for solving Lindblad MEs by using a Monte 
Carlo algorithm, even for multi level systems. Second, 
it is the only ME to guarantee complete positivity [30] , 
A crucial step to obtain a Lindblad ME is the secular 
approximation which neglects fast rotating terms. How- 
ever, its validity for time dependent systems is not gener- 
ally given and some work [2JH2E] uses non-secular MEs, 
which indeed result in negative probabilities for certain 
initial states. 

The aim of this work is to derive a completely positive 
ME for adiabatic, but otherwise very general systems 
(non-cyclic, multi-level), which is valid to exponential 
order e~ x l A in the adiabatic parameter. In particular, 
we rigorously justify the secular approximation leading to 
the desired Lindblad form. All other restrictions (weak 
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coupling, short bath correlation time, separable initial 
state) are identical to the time independent case [29] . 
Contrary to the time independent case, the Lindblad op- 
erators do not describe transitions between energy eigen- 
states, but transitions between super adiabatic states to 
be introduced below. The difference to the eigenstates 
is of order A. While in some instances this modification 
results in rather small quantitative effects and might be 
neglected (as in [2"2"H2~i] ). there are also situations (see 
our Landau-Zener example) where the super adiabatic 
states should be used, even for a correct qualitative un- 
derstanding. 

In adiabatic quantum computing, errors due to non- 
adiabaticities can be modelled by one or many succes- 
sive Landau-Zener (LZ) transitions. Therefore, we ap- 
ply our theory to the LZ problem where we find that 
pure dephasing does not alter the exponential decrease 
of the transition probability P g ^ e = e~ const / v with de- 
creasing transition velocity v, while the traditional Lind- 
blad ME 22 2.'? wrongly predicts transitions of order v 2 . 
Numerical solutions of our ME show that coupling to a 
low temperature environment even lowers the transition 
probability and stabilizes the adiabatic theorem. 

We now derive the central result of this paper: A mas- 
ter equation of Lindblad form governing the evolution 
of the density operator of an adiabatically driven sys- 
tem weakly coupled to a Markovian bath. We take the 
Hamiltonian of system and bath of the form 



Furthermore, we introduce 



Htot - H(t) + H h + A <g> B hl 



(1) 



where the subscript b refer to bath operators and no sub- 
script refers to system operators. The generalization to 
a general interaction along the lines of [29] is straight 
forward. Assuming a factorizing initial density opera- 
tor ptot(O) = p(0) (g) pb and a thermal state of the bath 
Pb oc e~ H b/ k BT ^ ^iq systems density operator in the in- 
teraction picture satisfies the master equation [29 



Pit) 



drTr[B(r)B Pb ] 



x [A{t-T)p(t)A(t)-A(t)A(t-T)p(t)]+h.c. (2) 



Here, the Born approximation ptotit) ~ pit) <S> ph and 
the Markov approximation p(t — r) w p{t) have been 
employed which are valid for weak coupling to a bath 
with short memory time Tb] and h.c. denotes the her- 
mitian conjugate. Operators with an argument are in- 
teraction picture operators, i.e. A(t) = W(t, 0)AU(t, 0) 
where U(t 2 ,ti) — Texp[— i J^dt'H(t')] is the time evo- 
lution operator of the closed system and T denotes time 
ordering. 

At each time we introduce a complete set of basis states 
Ifa), which we will later choose such that decoherence 
acts in this basis, i.e. the environment causes the den- 
sity operator to approach a diagonal state in this basis. 
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tft(t,t-r)|^)(3) 
(4) 



which arc, up to a phase factor, the states |<^) evolved 
backwards in time with |<^(i)) = \4>a)- The phase factor 
in Eq. ([3| is chosen such that the states are parallel trans- 
ported {4> t a {t~ T )\ £p l^ait — t)) = to minimize their 
rate of change with r. We neglect the change of the 



ifi r dt>Ei(t>) 



energies within and use exp 
Although our main results are valid without this approxi- 
mation, it is not a strong restriction for adiabatic systems 
and significantly simplifies the notations. 

Inserting twice the identity operator J2 a \(j>a) ((t>a\ for 



A(t) as well as EJ^i 
get 
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X {^{t-T)\A\4>%{t-T)){^ a ,\A\^ a ) 

x [C/t ( t , 0) 1 4) | U(t, 0)p(t)tf (t, 0) | fa )(fa\ U(t, 0) 

- S a prf (t, 0) |#,,)($, | U(t,0)p(t)] +h.c, (5) 

where we introduced the transition frequencies 0Jap, = 
{El, -El). 

At this stage we ask the question: Into which state 
does the system relax if the bath has zero temperature, 
i.e. if Ty[B{t)Bpb] has only positive frequencies? We 
can find a definite answer if (^(t — r)|v4|0L(i — t)) has 



only frequencies in r which are small compared to \uj 
Then, the r- integral of Eq. ( 5|) vanishes for Ep 



transitions /?<->/}' only go 



and 

owards lower energy. There- 



fore, the system evolves towards the state U^(t, O)|0^) 
(or \4>p) in the Schrodinger picture) with lowest energy, 
similar to an open quantum system with a time indepen- 
dent Hamiltonian. It is important to realize that this 
conclusion holds only if we can find a basis such 
that (<j)U [t-T)\A\^,(t-T)) with \fa(t-T)) as defined 
in Eq. Q does not oscillate with fast frequencies. 

For a time independent Hamiltonian one can trivially 
find a basis to satisfy the above statement, namely the 
eigenstates of the Hamiltonian. For adiabatic systems we 
can also find a basis which upon application of the ex- 
act time evolution operator changes only with frequen- 
cies much smaller than the transition frequencies, the 
so-called super adiabatic basis described below. 

If the system at some time t is in the instanta- 
neous eigenstate \n a (t)), then the time evolution oper- 
ator W(t, t — r) will cause the state to oscillate with the 
transition frequencies and with an amplitude of order of 
the adiabatic parameter A{t) defined by 



A 



max 

a/3 



(n a {t)\n p {t))\ 
\e a {t)-ep{t)\ 



(6) 
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FIG. 1: (Color online) The red line is the Bloch vector of the 
instantaneous ground state of the LZ problem and the dots 
are equidistant time steps. The Black curve is the third or- 
der super adiabatic ground state and the blue curve is the 
actual time evolution of an initial ground state. The time 
evolution operator causes oscillations around the super adi- 
abatic ground state. One result of this paper is that a low 
temperature environment causes the system to relax onto the 
black curve, rather than the red one, therefore stabilizing the 
adiabatic theorem. 



where e a (t) are the instantaneous eigenenergies of H(t) 
and the dot denotes the time derivative. If the sys- 
tem is in an eigenstate |n4(t)) of the first order super 
adiabatic Hamiltonian Hi = H + iRiR.\, then the os- 
cillation amplitude is only of order A 2 . The operator 
i?i = J2 a \ n a(0))(n a (t)\ transforms the system into the 
instantaneous basis. According to [I], one can use the 
j-th order super adiabatic basis \n 3 a (t)^ to further lower 
the oscillation amplitude to order A 3+ . However, once 
j sa 1/A the oscillation amplitude starts to diverge as jl; 
with the minimal amplitude of order e -1 ^ if the Hamil- 
tonian is analytic in time. An example is shown in Fig. [T] 

Therefore, in the following we make the explicit choice 
I'Pa) = \ n a{t)) where j is the natural number nearest to 
1/A and for adiabatic systems we can safely neglect the 
exponentially small oscillations caused by the evolution 
operator [31], [32]. The system will relax not into the 
instantaneous ground state, but into the super adiabatic 
ground state! We will see later that for LZ transitions 
this result is of importance for certain regimes. 

It is well known that the master equations ^ and ([5| 
do not preserve positivity |29) . The most common way to 
fix this is to perform a secular approximation which drops 
the terms in the interaction picture ME which rotate fast 
in t, as they average to zero. (Another method (33] uses 
a t- integral only from zero to t). The secular approxi- 
mation might give problems for time-dependent Hamil- 
tonians as many more frequencies become involved, some 
of them not large enough to justify the secular approx- 



imation. However, according to [3] H] W (t , 0) — 

|</>°)exp i f*dt' + O (e -1 /- 4 ) and upon neglecting 

the exponentially small non-adiabaticities, we are only 
left with the frequencies E l a = (0^| H(t) \4> a ), much like 
in the time independent case. We conclude that the 
secular approximation is justified for adiabatic systems 
if performed in the super adiabatic basis. Our final 
Schrodinger picture master equation of Lindblad form 
reads [M] 

p = -i[H + H hS ,p\+L oP Ll - ^{ijicp} 

+Y1 L opp L U - \{ L ip L a p,p} (7) 

with time dependent Lindblad operators 

Lo(t) = V^£(<^K)K)(^|, ( g ) 

a 

Lap® = \Jl{w a p) (4>a\ A \<f>p) |</>q)(0^| • (9) 

and the Lamb shift Hamiltonian H^s(t) = 

E Q ^(^)K^MI^>| 2 |^X^|- Here, 7 (w) and 
S(lu) are the real and imaginary parts of the Laplace 
transform of the bath correlation function Tr[B(T)Bpb], 
respectively. We note that neither Lamb shift i?LSi 
nor dephasing Lq affects the system if it is in a super 
adiabatic state, and that relaxations and excitations L a p 
change the populations of the super adiabatic states. 

A note on the time scales: Like all secular weak cou- 
pling master equations, Eq. Q - Q are valid provided 
the coupling is weak compared to the inverse bath corre- 
lation time t 6 -1 and that relaxation and dephasing times 
are long compared to the energy differences e a (t) — ep(t). 
Adiabaticity requires \(n a (t)\hp(t)) | <C \s a (t)— £p(t)\ and 
we also used \(n a (t)\hp(t))\ <C t^ 1 . We emphasize that 
we do not assume t^ 1 >• \e a (t) — ep(t)\, which would be 
valid only for extremely fast vanishing bath correlations. 

We now go on to apply our master equation Eq. Q to 
the Landau-Zener problem with 



H(t) 



-vt A 
A vt 



(10) 



Without coupling to a bath, the transition probability 
from the ground state at t = — oo to the excited state 
at t = oo was found [35H38] to be P g ^ e = e^ 2 /^). 
Although the sweep from — oo to oo is not realistic, finite 
sweeps through an avoided level crossing which are per- 
formed in many experiments are very well approximated 
by this transition formula. 

The coupling to the environment is via A — a Zl such 
that relaxation processes vanish in the limit t — > ±oo. 
In the case of a bosonic bath, there are very nice ex- 
act [521 HO] and numerically exact [JT] results, which 
show that a a z coupling to a zero temperature bath 
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FIG. 2: (Color online) The transition probability P 9 -» e for 
A = 1 as a function of 1/v without coupling (black), and 
with dephasive coupling 7(0) = 0.003 (red), 0.01 (green), 
0.03 (blue), and 0.1 (orange). In panel (a) we used Eq. |7| - 
|9| and the forth order super adiabatic states for |<^„), while 
in panel (b) we used the instantaneous eigenstates of H(t). 
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FIG. 3: (Color online) The transition probability P 9 ->e for 
A = 1 as a function of 1/v for an Ohmic environment. 
Panel (a) is for a low temperature environment T< A with 
70 = (black), 0.01 (red), and 0.03 (green). For panel (b) we 
used T = 0.5 and 70 = (black), 0.01 (red), and 0.1 (green). 



does not modify P g ^ e - The reason is an astonishing 
cancellation of amplitudes which move into excited bath 
states at different times. This strictly non-Markovian 
effect can not be captured by a master equation tech- 
nique applied here or in [21]. However, this cancellation 
of amplitudes requires that the temperature of the bath 
is exactly zero in a sense that all bath modes are in the 
ground state. While it is realistic in many experimen- 
tal setups to achieve temperatures small compared to 
the level splitting of the system, a temperature below 
the slowest bath modes is much more challenging [12] . 
Therefore, a Markovian weak coupling calculation might 
reveal interesting physics not captured by exact calcula- 
tions, as long as the bath correlation time is short (which 
is not the case for exactly zero temperature). 

In many solid state experiments the Hamiltonian can 
only be obtained qualitatively from first principles, while 
the precise level splittings are obtained experimentally 
and already include the Lamb-shift. For this reason we 
will discard i?LSj although for some applications involv- 
ing the use of a LZ transition to characterize the envi- 
ronment [24] . the Lamb-shift is of great importance. 

We obtained the transition probability P g -^ e by nu- 
merically solving Eq. ([7]), and plots versus the inverse 
transition velocity w _1 are shown in Fig. [2] (a) for differ- 
ent coupling strength to a slow environment, i.e. for pure 
dephasive coupling. The fact that all lines fall on top of 
each other show that dephasing has no effect on LZ tran- 
sitions. That is to be expected as a constant (or slow) 
environment coupled via a z merely shifts the time of the 
LZ transition. If we used the traditional ME 22 2.3 with 
Lindblad operators given in terms of instantaneous eigen- 
states, then we would not find an exponentially small 
Pg_>. e , but only a quadratically small one, as shown in 
Pig.|2j(b). This can be understood by using quantum 
trajectory theory |29j . A single trajectory continuously 
evolves along the blue line of Fig. [TJuntil a quantum jump 
rotates the state around the instantaneous ground state 
(red line). Subsequent continuous Hamiltonian evolu- 
tion leads to oscillations around the super adiabatic state 



(black curve). Averaging over different jump times lead 
to a population of the excited state of order of A~ 2 . The 
same argument also explains that pure dephasing does 
not alter the exponentially small transition amplitudes if 
the correct Lindblad operators are used. 

To discuss relaxation and excitation processes, we 
choose an ohmic bath correlation function 



with a cut-off frequency ui c = 5 A. Fig. [3] (a) shows the 
transition probability for a low temperature bath for var- 
ious coupling strengths 70. It can be seen that increas- 
ing the coupling is equivalent to decreasing the velocity 
v. The result that a low temperature bath can stabi- 
lize adiabatic ground state evolution is already known 
from Cooper pair pumping [25] , but contradicts the exact 
spin-boson calculations [SHI SO] ■ This is a manifestation 
of the reasoning described three paragraphs earlier. 

The finite temperature transition probability is plotted 
in Fig. [3](b). First, it decreases with 1/v due to improv- 
ing adiabaticity, and then increases because the environ- 
ment has more time to excite the system. For even slower 
velocities, in particular for stronger coupling, P g ^ e de- 
creases again, because the system has more time to relax 
at times when vt T . 

To summarize, we used the same premises which are 
generally used to derive a master equation of Lindblad 
form for time independent Hamiltonians, and extended 
the analyses to adiabatically time dependent Hamiltoni- 
ans. Our treatment shows that the validity of the secular 
approximation can be extended to adiabatic systems. We 
further found that the system decoheres into the super 
adiabatic states in the sense that off-diagonals of the 
density matrix diminish in this basis due to dephasing, 
and that the system tends to relax to the super adiabatic 
ground state if the environment has low temperature. We 
applied our ME to the Landau-Zcner problem where we 
found intuitive results while we also showed that a more 
naive ME with decoherence in the instantaneous eigen- 
states of the Hamiltonian would lead to wrong transition 
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probabilities. A positive result is that coupling to a low 
temperature bath can stabilize adiabaticity. 

Note: During the writing of this letter, a derivation of 
a quite general adiabatic ME with secular approximation 
appeared [33] , which describes approximations similar to 
here in much detail and is therefore of great value. The 
main difference to our work is that they assume U'(t, i — 
t) s» e lT - H (*) ! similar to [24], and therefore do not use 
super adiabatic states in the Lindblad operators. As is 
clear from their and our work, this neglects contributions 
of order ^A (unless t^ 1 >• w a p)- 
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